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This paper discusses the analytical solution of a finite circular plate 
on an elastic foundation centrally loaded by the rigid indenter; and the 
procedure to use NASTRAN as a subroutine to iteratively converge to this 
solution numerically. 


INTRODUCTION 


In contact recording applications, where a thin flexible spinning disk 
backed by an elastic pad is penetrated by a read/write head, it is 
important to understand the head penetration, the contact area, and the 
pressure distribution. A computer program was written to iteratively 
converge to the solution using MSC/NASTRAN version 60* (hitherto called 
NASTRAN) as a subroutine. Results were erratic. It was found that 
extremely fine grid was needed around the contact area to converge to 
the right solution. 

In order to calculate the approximate contact area, the model was 
simplified to a finite circular plate on an elastic foundation centrally 
loaded by the head. The approximation is fairly valid because the 
flexural rigidity of the plate is extremely low and the reaction forces 
at the center of the spinning disk is about 2% of the total load. The 
effect of most of the load, thus, is to cause local deformation of the 
plate between the head and the pad. An analytical solution for this 
case was obtained. The iterative procedure to converge to this solution 
by using NASTRAN as a subroutine to verify the analytical solution was 
also defined. 


MacNeal-Schwendler Corp version of NAsa STRuctural ANalysis. 
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ANALYTICAL SOLUTION 


Consider in figure 1 the spherical indenter with radius being pressed 

into the circular plate with center at the origin, the plate resting on 
an elastic foundation with elastic spring constant k. It is of interest 
to determine the inside contact radius OA and the pressure distribution 
on the indenter over that contact. Notice, as shown in figure 1, that 
the outside contact radius is OB and the plate physically separates from 
the elastic foundation for radii larger than OB. Behavior of this 
"free" plate can also be modeled easily by matching boundary conditions 
at B, but it was not necessary in our problem. Instead we were seeking 
the solution to the problem "what is the outside radius R q of the plate 

such that vertical displacement W=0 at R^?" 

The behavior of the plate on elastic foundation can be described (ref. 1) 
by equation 
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For our case, in which the indenter presses on the center of the plate, 
there is a radial symmetry, and the solution is independent of angle 0, 
we can rewrite the above equation as: 
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Where W = vertical displacement 

q = distributed lateral load 
D = flexural rigidity of the plate. 
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h = thickness of the plate 
E = Young’s modules of the plate 
v = Poisson' s ratio of the plate 

k = spring rate/unit area of the elastic foundation. 
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In seeking the solution for above, we divide the plate into two sections, 
as shown in figure 2. The inside (section 2) portion conforms to the 
sphere up to contact radius R^. The outside portion of the plate must 

also satisfy differential equation (1) with appropriate matched boundary 
conditions . 


OUTSIDE PORTION OF THE DISK 


Consider now the outside portion of the disk (section 1 in figure 2) . 

If q = constant, such as gravity load, the particular solution of 
equation (1) is simply Therefore, we now seek the complementary 
solution (in absence of q) of 
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where R is the radius of the indenter. 
s 

Equation 2 can then be rewritten as 
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This is a linear differential equation of the 4th order. The solution 
can be represented in terms of Bessel's function. However, it is 
presented below in terms of series solution for ease of computation. 
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The general solution is given by 


n = A 1 «1 + A 2 «2 + A 3 <$ 3 + A 4 6 4 
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where A. . . A, are constant, 
1 4 
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Coefficients in above two series can be generated by recursive relationship 
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where 
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A = _J_ C 6 1540 C 10 (7) 
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and coefficients in 6 and 7 can be recursively generated by 
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Four constant A^..A^ can be thus evaluated by following boundary conditions. 
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iortion of the disk can be obtained 
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SOLUTION OF INSIDE PORTION OF DISK 


The solution for the inside portion of the disk (section 2, figure 2) 
that conforms to the sphere satisfies 
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where VL = deflection of the disk at r = R^ obtained from equation (3) 


PRESSURE DISTRIBUTION 


Total force F is given by 
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Pressure distribution under the head is given by equation (8) multiplied 
by k. However, at radius R CI there is additional shear force required to 

keep the outside portion of the disk in equilibrium, and is given by 
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NORMALIZING 


Solution to the equations are functions of R g , F, £, R q . However, if 

the disk is very weak and outside radius R q is taken to be that value 

for which W=0, then from equations (3 through 8), it can be shown that 

FR , R A , OR £ , qR £^ are all normalized functions of 

s * o' x s ’ ■ s 

R ci/r 

Figures 3 and 4 show some of the relationships. 

NASTRAN MODEL DESCRIPTION 


The application of the static analysis of MSC/NASTRAN to the problem in 
figure 1 is demonstrated in this section. 

The elastic foundation in figure 1 is first discretized into the "gap 
scalar spring" of finite length corresponding to grid points on the 
plate. The "gap scalar spring" is defined as a linear spring that can 
be compressed only. An extension of this spring will produce no spring 
force. This spring force will then be used to generate the FORCE card 
for the NASTRAN program. Secondly, an isoparametric bending element 
with transverse shear is used to model the plate between the indenter 
and the elastic foundation. As shown in figure 5, to save computer 
time, only a section of this circular plate was modeled in this work. A 
symmetric condition is used along the edge boundaries. Since the spherical 
indenter can be treated as a rigid body in this problem, its contour is 
computed and stored in [C] . This [C] is used as an enforced displacement 
to the plate model within the contact area. 

The radius R of the contact area between the indenter and the plate is 

L* X 

essentially the solution we need. In the following part, an iterative 
technique was developed to compute the R^, contact radius, and the load 

distribution over the plate. 

Final solution to the problem must meet the following conditions: 

1 There should be no geometric interference between the deformed 

plate and the contour of the indenter. The interference example is 
as shown in figure 6, where the plate deformation pattern IW]^ 

interferes with the indenter contour [C] . The smallest R^ without 

a geometric interference is the solution to this problem. 
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2 The s umma tion of the distributed load, or the "gap scalar spring" 
force, over the plate be equal to the total force F. 

£ [q] [area] = F (9) 

In order to satisfy the above two conditions, an iterative type of 
algorithm was developed to solve the problem by utilizing the CALL 
NASTRAN technique. In this case, the whole NASTRAN program becomes a 
subroutine that can be implemented by any user-developed program. The 
NASTRAN program can thus be invoked by using the CALL statement in the 
user's program just like using any other conventional subroutines. 

To start this iteration process, an initial guess of R CI is required as 

an input to the program. With this assumed R , the static contact area 

between the indenter and the plate can be defined and set equal to the 
indenter contour by using the enforced displacement card, SPC. Initial 
plate deflections are obtained by using NASTRAN. With these deflections 
and the elastic foundation stiffness k, forces on the grid can be calculated 
satisfying condition (2) above. This [q], as the FORCE cards, along 
with the previously defined SPC cards, are then put into the NASTRAN 
program to compute the plate displacement [W] . This plate displacement 
is not supposed to interfere with the indenter contour IC] . If this 
geometric constraint is not satisfied, a larger R CI will be assumed next 

and repeat the above steps until a R^ can be found that satisfies this 

geometric constraint. This final R^ is the solution we are trying to 

obtain. Thus, the static contact area and the final distributed load 
can be computed accordingly without any difficulty. 

Figure 7 shows the simplified flow diagram of this iteration process. 

Some additional explanations and programming considerations are made 
below. 

1 Since the NASTRAN program is being called and invoked many times, 
NASTRAN CHECKPOINT and RESTART features have been used. 

2 The CALL NASTRAN statement actually includes the following functions: 

a Generate the NASTRAN executive control deck, including CHECKPOINT 
dictionary for RESTART use. Case Control Decks and the Bulk 
Data Card. 

b Invoke the NASTRAN program by using a CALL statement. The 
NASTRAN program is invoked throughout the entry point of its 
load module. This entry point name is "NASTRAN" in our case. 
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c Recover the displacement vector computed by NASTRAN as the 
output from this GALL NASTRAN operation. The module OUTPUT4 
was used to modify the DMAP sequence to get the desired output 
from NASTRAN. 


RESULTS 


The iteration scheme defined in the NASTRAN model gave results quite 

close to those obtained analytically. The graphs are shown in figures 3 

and 4. Given force F and the radius R of the indenter, contact radius 

s 

and the pressure distribution can be obtained from figures 3 and 4 
respectively. 
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SPHERICAL INDENTER WITH PLATE ON ELASTIC FOUNDATION 


FORCE = F 



Figure 1. Spherical indenter with plate on elastic foundation 



PLATE SHOWN IN TWO SECTIONS. INSIDE PORTION O-A, 
CONFORMS TO THE SPHERE WITH RADIUS Rg. 


SECTION 1 



Figure 2. Plate shown in two sections. Inside portion, O-A, 
conforms to the sphere with radius R s . 




Figure 3. Contact radius 











SECTION OF THE CIRCULAR PLATE MODEL 



i Figure 5. A section of the circular plate model 
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SIMPLIFIED FLOW DIAGRAM OF THE ITERATION PROCESS 



'Figure 7. Simplified flow diagram of the iteration process 






